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Abstract. This paper deals with higher gradient integrability for cr-harmonic func- 
tions u with discontinuous coefficients a, i.e. weak solutions of div(trVu) = 0. We 
focus on two-phase conductivities a : !1 C I 2 4 {01,0-2} C M 2x2 , and study the 
higher integrability of the corresponding gradient field |Vu|. The gradient field and 
its integrability clearly depend on the geometry, i.e., on the phases arrangement de- 
scribed by the sets Ei = a^ 1 (ai). We find the optimal integrability exponent of the 
gradient field corresponding to any pair {ci, 02} of positive definite matrices, i.e., the 
worst among all possible microgeometries. We also show that it is attained by so- 
called exact solutions of the corresponding PDE. Furthermore, among all two-phase 
conductivities with fixed ellipticity, we characterize those that correspond to the worse 
integrability. 
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1. Introduction 

Let H be a bounded, open and simply connected subset of M 2 with Lipschitz continuous 
boundary. We are interested in elliptic equations in divergence form with L°° coefficients, 
specifically, 

(1.1) div(o-Vu) = in 

1 
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Here a is a matrix valued coefficient, referred to as conductivity, and any weak solution 
u G Hi oc (fl) to the equation is called a a -harmonic function. The case of discontinuous 
conductivities a is particularly relevant in the context of non homogeneous and composite 
materials. With this motivation, we only assume ellipticity. Denote by M 2x2 the space of 
real 2x2 matrices and by M 2 *,^ the subspace of symmetric matrices. 

Definition 1.1. Let A G (0, 1]. We say that a G L°°(ft;M 2x2 ) belongs to the class M{X,Q) 
if it satisfies the following uniform bounds 

(1.2) cr£-£>A|£| 2 for every £ G M 2 and for a.e. x G n, 

(1.3) cr _1 £ • £ > A|£| 2 for every £ eR 2 and for a.e. x £ Q, 
Moreover, we denote by M sym (A, f2) i/ie sei M 2 ^ n M (A, O). 

The reader may wonder why to use the notion of ellipticity given in Definition 11.11 For an 
explanation related to its relationship with if-convergence we refer the reader to [2]. 

It is well known that the gradient of cr-harmonic functions locally belongs to some L p 
with p > 2. The main goal of this paper is to explore this issue, focusing on two-phase 
conductivities {ci, c 2 } C M. We will review known results and prove some new one. 

Any cr-harmonic function u can be seen as the real part of a complex map / : n \— > C which 
is a H\ oc solution to the Beltrami equation 

(1.4) fz = nfz + vTz, infi, 

where the so called complex dilatations \l and v, both belonging to L°°(n,C), are given by 

o- 22 - an - i((?l2 + 0"2i) 1 - det cj + i(cri2 - cr 2 i) 

(1-5) u = , v = , 

v ; p l + tr<T + detc7 l + trcr + deter 

and satisfy the ellipticity condition 

(1.6) ||H + HH^oo < 1. 

Let us recall that weak solutions to (|1.4p are called quasiregular mappings. They are called 
quasiconformal if, in addition, they are injective. The ellipticity (|1.6p can be expressed by 

K-l 



(1-7) II H + HIU- < J , - 

for some K > 1. The corresponding solutions to (jl.4p are called K-quasiregular, and X- 
quasiconformal if, in addition, they are injective. In 1994, K. Astala [3] proved one of the 
most important pending conjectures in the field, namely that planar K- quasiregular mappings 
have Jacobian determinant in l \ Astala's work represented a benchmark for the issue 

of determining the optimal integrability exponent which was previously studied in the work 
of Bojarski [7] and N. Meyers [13] . 

Summarizing, to any given cr G M (A, n) one can associate a corresponding pair of complex 
dilations via (|1.5p and therefore, via the Beltrami equation (|1.4p a quasiregular mapping. 
Therefore, given A G [0, 1) and given a G M. (A, 0) one can find K = K(a) by using (11. 5p and 
(11. 7p in such a way that the cr-harmonic function u, solution to (jl.ip is the real part of a K- 
quasiregular mapping. The Astala regularity results in this context reads as |Vu| G L^ eak [Q), 
where p K ■= 

A more refined issue is to determine weighted estimates for the Jacobian determinant of 
a quasiconformal mapping. A first result in this direction was given in [6]. A much finer 
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recent result, is given in [3], see formula (1.6). Throughout the present paper we focus on the 
simpler framework of L p spaces. 

The first question is to determine the best possible (i.e. the minimal) constant K{a) such 
that if u is a-harmonic with a G Ai (A, fi), then u is the real part of a i^(cr)-quasiregular 
mapping. Astala writes in his celebrated paper that his result implies sharp exponents of 
integrability for the gradient of solutions of planar elliptic pdes of the form (jl.ip . and he says: 
"note that the dilation of / and so necessarily the optimal integrability exponent depends in 
a complicated manner on all the entries of the matrix a rather than just on its ellipticity" . 
Alessandrini and Nesi [2], in the process of proving the G-stability of Beltrami equations, 
made a progress which can be found in their Proposition 1.8. Let us rephrase it here. See 
also [1] for the estimate fjl .9|) . 

Proposition 1.2. Let A G (0,1]. Then 



(1.8) K X := sup K{a) = I±^_L_i! 

aeM(X,Q) A 

(1.9) K s x ym := sup K{a) 



1 

A 



In Section [2]2] we give a simpler and more geometrical proof of Proposition II . 21 based on the 
real formulation of the Beltrami equation (see Propositions 12.21 and 12.3]) . In [2], pp. 63, the 
authors noticed that, the supremum in (|1.8p is attained on specific non symmetric matrices. 
As a straightforward corollary, in [2j the authors write the version of Astala's theorem which 
is adequate for matrices belonging to Ai (A, Q) that we recall here in an informal way. Any 
o"-harmonic function with a G A4(X,fl) satisfies the property |Vit| G L K, where K\ is 
given by (jl.8p and px x '■= k^-i • This has to be compared with the version that holds true 
assuming a priori that a G Ai sym (A, O). In that case K\ can be replaced by K sym defined in 
(jl.9p . Optimality in the latter case was proved by Leonetti and Nesi |12| which began their 
work using the bound (jl.9p which had been already observed in Alessandrini and Magnanini 
[I]. Optimality means that there exists c G J\A S y m (\, fl) for which the estimate | Vit| G L w J^p, 
is sharp. 

Later there has been a number of increasingly refined results showing optimality of Astala's 
theorem for a different class of symmetric matrices a. Specifically Faraco [8] treats the case 
of two isotropic materials, i.e. when a takes values only in the set of two matrices of the form 
{KI, with I the identity matrix, which was originally conjectured to be optimal for the 
exponent -j^Sf by Milton [15]. In a further advance a more refined version was given in [5], 
where the authors proved optimality in the stronger sense of exact solutions. 

However the original question implicitly raised by Astala was apparently forgotten. In 
this paper we go back to that and we prove optimality for a generic two-phase matrix field 
a G Ai(X,Q). To describe our approach let us first recall that when a is smooth, the 
corresponding cr-harmonic function is necessarily smooth and hence with bounded gradient. 
So the issue of higher exponent of integrability is really related to discontinuous coefficients. 
The simplest class of examples is when one has a conductivity taking only two values. We 
therefore ask the following questions. Given two positive definite matrices, o\ and o"2, consider 
the class of matrices a G At (A, Q) of the special form a(x) = o\xe y +&2Xe 2 i where {E\, E2} 
is a measurable partition of f2 and XE t denotes the characteristic function of the set E{. In 
the jargon of composite materials this is called a two-phase composite. What is the best 
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possible information one can extrapolate from Astala's Theorem? As already explained, to 
the ellipticity A of a there corresponds a suitable constant K (cr) in the Beltrami equation. We 
are naturally led to the following related question: given //, v £ L°°(£l; C) satisfying (|1.7p with 
K(/j,, v) > 1, is it possible to transform /i and v, by a suitable change of variables, specifically, 
by affine transformations, in order to decrease K and thus gain a better integr ability for 
the solution of the transformed Beltrami equation? The key observation here is that the 
summability of solutions of the Beltrami equation is invariant under such transformations, 
while K(fj,,f) is not. It is then well defined the minimal Beltrami constant K' mm attainable 
under such transformations. In Proposition 15.41 we find an explicit formula for such K mm 
in terms of all the entries of o\ and oi- Moreover, K mm gives a sharp measure of the 
integrability properties of solutions to (jl.ip . This is stated in Theorem 15.11 which, for the 
reader's convenience, we reformulate here in a more informal way. 

Theorem 1.3. i) Let a G .M(A,S1) with a G {ci,C2}- Every a-harmonic function u satisfies 
Vu £ Lf oc (£l) for every p £ [2,p Kmi n). ii) There exist a G with a £ {ci,^} and a 

a-harmonic function with affine boundary conditions such that, for every ball B C f2 

(1.10) / \Vu\ p K™ n dx = oo . 

Jb 

A key step to prove Theorem 11.31 is to prove the optimality of Astala's Theorem for a new 
class of symmetric conductivities, specifically, for matrices of the form 

(1.11) ( r = XE 1 diag(S 1 ,A- 1 ) + Xi? 2 diag(S 2 ,A), with A < Si, 5 2 < A -1 , 

thus generalizing the isotropic case Si = A , S2 = A, considered in [5] and [3]. 

As a corollary of Theorem II. 3} we prove that the bound (|1.8j) for non symmetric matrices 
too is optimal. Indeed, there exists a £ A4(X, tt) of the form 

(1.12) <* = XE 1 (" b b a )+XE 2 (l witha = A, b = ±VT^V 

and a c-harmonic function u such that the bound Vu £ L \ is sharp (see Theorem 15. 2[) . 

Finally, a natural question, both in the symmetric and in the non symmetric case, is wether 
there are other two-phase critical coefficients, that is to say, two-phase coefficients a for which 
the bounds in Proposition II .21 are attained and optimal in the sense of (|1.10p . In Theorem 15. 31 
we give a complete answer to this question, characterizing all the critical conductivities with 
fixed ellipticity. In the symmetric case, the critical conductivities are given (up to rotations) 
exactly by those in (jl.lip (for suitable partitions E\, E<i). In the non symmetric case, the 
only critical conductivities are as in (|1.12p . 

We remark that one can find optimal microgeometries for cr's which are not two-phase. 
The simplest example is given by a "polycrystal" like in the first example given in Leonetti 
and Nesi [12]. In that case a is symmetric, the eigenvalues are A and A -1 but the eigenvectors 
change from point to point. 

2. More about ct-harmonic functions and the Beltrami system 

In the present section we review some well-known connections between cr-harmonic functions 
and the Beltrami system which we use in the rest of the paper. We refer the interested reader 
to [2] for a more detailed presentation of the argument. 
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2.1. Complex vs real formulation of a Beltrami system. Consider the Beltrami equa- 
tion (jl.4p . It can be rewritten in the equivalent form 

(2.1) Df l HDf = G det Df , 

where G and H are real matrix fields depending on \i and v. Specifically, 

1 / |l + /u| 2 - \v\ 2 23f(/x) 



(2 ' 2) °~d\ 23( M ) |1- M | 2 -|^| 2 

1 / |1_ H 2_ H 2 _ 2QM 

d V -29(f) |l + i/| 2 "H 2 

where 



d = V(i-(H-H) 2 )(i-(H + H) 2 )- 

We will refer to (|1 .4[) as well as to ()2.ip as the Beltrami system. Let SL{2) be the subset of 
M 2x2 of the invertible matrices with determinant one, and let SL sym (2) = M 2 *^ n SL(2). 
Notice that G and H belong to SL sym (2) and they are positive definite. In fact injective 
solutions to ()2.ip have a very neat geometrical interpretation. They are mapping / : O, — >■ O' 
which are conformal, i.e., they preserves angles, provided one uses the right scalar products, 
namely the one induced by G in £1 and H in f2'. This interpretation has many consequences. 
We will get back to this point later in the paper. Inversion of the above formulas yields 

_ Gn — G22 + 2iG\2 _ H22 — Hn — 2iH\2 

G\\ + G22 + + -£^22 G11 + G22 + .H11 + -£^22 

By combining (|2.2p and (|1.5p we obtain a formula for G and H as functions of a, 
(2.3) 

1 /- a22 _£i2±£2i \ tw % 1 /- det a ~ CT12 +' 721 



G(a) " TleTP I J ' H{a) ~ TdetP I 1 



where a 5 = . Inversion of (|2.3p gives 



(2.4) < 7=-J-(G- 1 + fTi 2 j) 



H22 

' -1 



where 

( 2 - 5 ) J 1 1 „ 

Moreover, we can express <r as a function of /i, f inverting the algebraic system (11.51) . 



/ |1-/.| 2 -|H 2 2Q^- M ) 
/ IH 



(2.6) a 



V 



|l +i ,p_| M |2 | 1+1/ |a_| M |2 
-2S(i/+/z) jl+^-H 2 

|i+H 2 -IH 2 |i+H 2 -IH 2 



Let us clarify the relationship between the Beltrami equation and cr-harmonic maps. Given 
positive definite matrices G and H in L°°(Q; SL sym (2)), let / = (u, v) be solution to (|2.ip . 
Then, the function u is cr-harmonic, with a defined by (I2.4p . Conversely, given cr satisfying 
the ellipticity conditions (jl.2p - (jl.3p and given a cr-harmonic function u, the map / := (u, v) 
solves (|2.ip . where G and are defined by (j2.3|) . u is such that 

(2.7) J T Vt> = crVu, 
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and J T is the transpose of J denned in (|2.5p . The function v is called stream function of u, 
and is defined up to additive constants. Moreover, ||V/||lp is finite if and only if ||Vu||lp is 
finite. 

2.2. Different formulations of ellipticity and higher gradient integrability. Here we 
introduce classical notions of ellipticity for elliptic and Beltrami equations, and we recall 
the fundamental summability results due to Astala [3] and some of its consequences due to 
Leonetti and Nesi |12j . From now on, we will always assume that the values of /u , v ,G , H 
and a are related according to (jl.5p and (|2.2|) . 

The ellipticity corresponding to any pair ^i, v G L°°(£l;C) satisfying (jl.6p is the positive 
constant k(fi,v) defined by 

(2.8) k([i,v) := + \v\\\ L °o. 

An alternative measure of ellipticity, that will be most convenient in our analysis, is provided 
by the following quantity 

Having in mind (|2.2p . we define k(G,H) and K(G,H) in the obvious way, i.e., 
(2.10) k{G,H) = k(ji,u), K(G,H)=K(n,v), 

and whenever no confusion may arise, we will omit the dependance on their argument. In the 
next proposition we give a more explicit formula for such ellipticity. We will denote by g(x) 
e h(x) the maximum eigenvalue of G{x) and H(x), respectively. 

Proposition 2.1. Let G, H G L°°(J7; SL sym {2)) be positive definite. Then 

(2.H) K =\\gh\\ L oo {n) . 

Proof. A direct computation shows that the maximum eigenvalues of G and H are given by 



^(i-H + H)(i + H + M) , v^ + M-HXi + M + H) 

9= rj- : , , =F , , , == , h~ 



V(i + H-|H))(i-H-IH)) vFR + WTl 

Therefore qh = — \ L , which yields 

i-CN + M) 

\\ 9 h\\ L ~ = 1+ ^\ +H }r = 1 -±4= K . 

ny l|L i- H + M x i-fe 



□ 



Next, we relate the ellipticity bounds for the second order elliptic operator (11. lj) with 
the ellipticity of the associated Beltrami equation. Following the notation of (I2.10p . we set 
K(a) := K(G,H), where G,H and a are related by (f2T3]) - ([2~4|> . The following result has 
been proved in [12] and [2]; for the reader's convenience, we give here a proof based on 
Proposition 12.11 



Proposition 2.2. Let A G (0, 1]. For each a G M. (A, fi) we /lave 
(2.12) A - M <1±4ZZ. 
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If in addition a is symmetric, then 

(2.13) K(a)<j. 

Proof. Let Ai,A2 be the eigenvalues of a s , with Ai < A2. Then, from the assumption a G 
Ai (A, ft) and the relationship 

, s deter 5 g l 

( CT ) = -7-7 — ( CT ) , 
det a 

it follows 

(2.14) A 2 >Ai>A, 

to 10 det a s Ai 

A2 det cr det a 

Next let g and h be the largest eigenvalue of G and respectively. By (|2.4p . it is readily 
seen that 



cr 



s 



1 



■G- 1 . 



and hence 
(2.16) 

From (1231) it follows 



#22 

1 1 Vdet a s 



H22 Ai Ai 



(2 - 17) ft+ s = CT (de " +1) ' 

Set P := ^ = - . Solving (|2.17p and choosing the root which is bigger than one, yields 
V det a s 



p _L ,/p2 _ 4 

(2.18) h= — V -. 

v 7 2 

Then, using (|2.16|) - f|2.18j) and the inequalities (|2.14p - ()2.15p . we obtain the following upper 
bound for gh 



9h= 2X[ 



det a + 1 + J (det a + l) 2 - 4 det er £ 



< 



2Ai 



T + 1 + VvT +1 



4A 2 



1/1 1 1 1 

77 V + V- + Wt77 + To + 



2 \ A A 



A 2 Af AAi 



A 
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Now suppose that a is symmetric and denote by Ai and A2 its eigenvalues, with Ai < A 2 . 
Since a G M. (A, t, , we have 



(2.19) 

Formula (|2.3p reduces itself to 

G = Vdetaa' 1 

Therefore 

(2.20) g 
In the case when A1A2 < 1, we find 



A < Ai < A 2 < T 
A 



H 



1 



— V det a, 
Ai 



If otherwise A1A2 > 1, we have 



K 



K 



1 



V det a 
1 

V det a 
1 

< - 

l°° A 



deter 
1 



max{AiA2, !}• 



IA2HL 00 < -r ■ 



□ 



In the next Proposition we look at conductivities a attaining the bounds ()2. 12j) and (|2,13p . 

Proposition 2.3. Let a £ Ai (A,f2) for some A E (0, 1]. Then the bound (|2,12p is attained 
if and only if on a set of positive measure there holds 

(2.21) a = 



a b 
—b a 



^ , with a = A, b = ±\/ 1 — A 2 



Moreover, if a is symmetric f)2. 13|) is attained if and only if either (| 1 . 2 j) or (jl.3p is attained 
on a set of positive measure. 



Proof. Keeping the notation introduced in the proof of Proposition 12.21 one can see that the 
bound ()2.12p is attained if and only if the inequalities (|2.14p - (|2.15p hold as equalities, namely, 

Ai 



A, 



Ai = A, 



det a 



A. 



It is readily seen that this is equivalent to (|2.2ip . The symmetric case is left to the reader. □ 

We now recall the higher integrability results for gradients of solutions to (11. ip and (|1.4jl . 
2K 

For K > 1, set px '■= — t- We start with the celebrated result in [3]. 
K — 1 

Theorem 2.4. Let f £ .H^QjC) be solution to ([HI) with K(n,v) > 1. Then 

V/eLf oc (0) V P €[2,p KM ). 

Recall that K\ and K^ m are defined by (jl.8p and (|1.9|) . respectively. A straightforward 
computation yields 



2 + 2^/T^A 2 " 



(2.22) PKs = , V K sym 

V 7 l-A + v 7 !^ ^ 1 + A 

As a consequence of Proposition 12.21 and Theorem \2A\ we obtain the following result which 
was proved in [12], [2]. 
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Theorem 2.5. Let a G A4(X,Q) for some A G (0,1]. TTjen, any solution u G f/j^ c (Q) to 
(II. ip satisfies 

VuGLf oc (Q) Vpe[2,p Kx ), 

and, ifaeL°°(n;M 2 s *£), 

Vu£Ll c (n) VpG [2,p K sy 
where px x and Px s x ym are given in (I2.22D . 

We are now ready to perform linear change of variables both in the domain and in the 
target space. It will be convenient to work with the real formulation of the equation (12. ip . 
Let A, B G SL(2) and set 

(2.23) f{x) := G(x) := B t G(Bx)B, H(x) := A t H(Bx)A. 

A straightforward computation shows that, whenever / : £1 \— > M 2 is solution to ()2.ip . / solves 

(2.24) DpHDf = G det Df. 

Clearly / enjoys the same integrability properties as /. This motivates the following definition, 

(2.25) K min (G,H) := min \\g(A, B)h(A, B)\\ L ~ , 

A,BgSL(2) 

where g(A, B) and h(A, B) denote the maximum eigenvalue of G and H, respectively. Remark 
that g(A) > c|| A|| 2 , h(B) > c\\B\\ 2 . Therefore, the minimum in (|2.25p is attained. Recalling 
()2.1ip . a straightforward generalization of Theorem 12.41 leads to the following result. 

Proposition 2.6. Let G,H G SL sym {2) and let K min (G,H) be defined as in (l2T25l . Then 
any f G H\ oc (Q; M 2 ) solution to (|2.ip satisfies 

V/GLf oc (0) VpG %p Kmin{G>H) ) . 

Remark 2.7. From the point of view of a -harmonic maps, Proposition ^. 6\ may be rephrased 
by saying that any solution u G Hy oc (Q) to (jl.ip . satisfies 

VueLfjn) V P e[2,p Kmin{a) ), 

where K mm (a) is defined in the obvious way, i.e., K m ' in {a) := K mm (G,H), and G,H and a 
are related by (j2.3|) . 



3. Examples of weak solutions with critical integrability properties 

In [8] , [5] , the authors exhibit an example of weak solution to (jl.ip with critical integrability 
properties. In their construction the essential range of a consists of only two isotropic matrices, 
namely, a : $7 i— > {K~ 1 L, KL} with K > 1. In this section we generalize their construction to 
the case 

(3.1) a 1 :=diag(K,S , 1 ), a 2 := diag^" 1 , 5 2 ), \z < S t < K, 

thus proving optimality of Astala's theorem for the whole class of matrices above. In Section 
Owe will show that such class cannot be further enlarged. 
We will need the following definition. 

Definition 3.1. The family of laminates of finite order is the smallest family of probability 
measures £(M 2x2 ) on M 2x2 such that 
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(i) £(M 2x2 ) contains all Dirac masses; 

(ii) if £™ =1 a i § Ai e £(M 2x2 ) and Ai = aB + (1 - a)C with rank(5 - C) = 1, then the 
probability measure ^™ =2 a^Ai + otiiptB + (1 — a)C) is also contained in £(M 2x2 ). 

Given v G £(M 2x2 ), we define the barycenter t> of v as 

P := / Mdv(M) . 

Theorem 3.2. Let o~\,o~2 be as in f)3. 1 1) . There exists a measurable matrix field a : ft — > 
{ci)C2} sitc/i t/iai i/ie solution u G io 

^ I div(crVu) =0 in Q 

[ u(x) = x\ on dft 

satisfies for every ball B C ft, 

/ \Vu\ PK dx = oo. 

The proof follows the strategy in [5\ Theorem 3.13], where the result is proved for o~\ = KI, 
o"2 = K~ l I. Here the main difference is that we work with coefficients that are not isotropic. 
For the reader's convenience we shortly reproduce the arguments of [5] pointing out the 
essential modifications. 

Step 1 {Reformulation of (|3.2[) as a differential inclusion). Recall that u is solution to 
(13. 2p if and only if u = /i where / = (/i, /a) is solution to the associated Beltrami equation. 
It is easily checked that, for a of the form (13. the latter condition is equivalent to 



(3.3) Df G E := E x U E 2 , 

where 

* = { {jIt) < t g r2 } ' ^ 2 = { (jJ 2 t) < t g r2 

The goal is to find a solution / G i^T 1 K. 2 ) to the differential inclusion (13. 3|) satisfying in 
addition the boundary condition f\{x) = x\ on dVt. 

Next we define a setting where to apply the Baire category method. Fix 5 > such that 



(«) * < f ( !' 1/ ^ )( 5r' ) ) 

V 4max{Si,5 2 }# 2 / 

and let 



(3.5) E:=En\( au au G M 2x2 : |a 12 | < 8a n 

Notice that the introduction of the small parameter 5 enforces the solutions to have gradi- 
ent pointing in a direction relatively close to /. This property hides the anisotropy of the 
coefficients <7j, and allows us to follow the strategy of [5]. Define IA as the interior of the 
quasiconvex hull of E (defined as the set of range of weak limits in L 2 of solutions to (|3.3p ). 
The following characterization of U holds 

E 10 ' 1 =U = E pc , 
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where E > and E pc denote the first lamination hull and the polyconvex hull of E, respectively. 
We refer to [H Lemma 3.5] for the proof of the identity above and for the notion of first 
lamination hull and polyconvex hull. 

Set _ 

X = {/ £ W lt0 °(Sl;R 2 ) : f piecewise affine,D/ 6 U a.e., /| 9n = x}, 
let X be its closure in the weak topology of H 1 , and denote by (X, w) the set X endowed 
with the weak topology w of H . Remark that I 6 U and therefore the set X is not empty 
as it contains the map f(x) = x. 

Step 2 (Existence of solutions by the Baire category method). The existence of solutions 
to the differential inclusion is proved by an application of the Baire category method, and is 
based on the fact that the gradient operator D : X i— > L 2 (S7;M 2x2 ) is a Baire-1 mapping, i.e., 
the pointwise limit of continuous mappings. We refer to [Tfll page 57] and references therein 
for further clarifications on this subject. The existence result is stated in the next theorem. 
We refer to [5l Lemma 3.7] for its proof. 

Theorem 3.3. The space (X, w) is compact and metrizable. Each f £ X satisfies f £ IA and 
f\&0, = x. The metric d on X is equivalent to the metrics induced by the L 2 and L°° norms. 
Moreover, the points of continuity of the map D : (X,w) — > L 2 (il; M 2x2 ) form a residual set 
in (X, w). Finally, any point of continuity f £ X of D satisfies Df £ E\ U E2. 

We deduce that the set of solutions to the differential inclusion f)3.3f) is residual in (X,w). 
The proof of Theorem l3.2l is then a consequence of Theorem 13.31 and of the following theorem. 

Theorem 3.4. The set 

is residual in X. 

Theorem 13.41 is proved following the same strategy of the proof of Corollary 3.12 in [Sj. We 
recall that in [5] the isotropic case S\ = K , S2 = 1/K is considered. In the present setting 
the proof is identical except for the proof of a key ingredient (namely, [5] Proposition 3.10]). 
Therefore, we only state and proof such result in Lemma 13.51 below. For this purpose it is 
convenient to introduce some notation. Given a matrix A = (a%j) £ M 2x2 , we denote by 
and A a its diagonal and anti-diagonal part, namely 

m * - ( 7 I ) • — ( I 7 ) • 

Moreover we will identify Ad and A a with points of M 2 : Ad = (011,022); A a = (012,021). 
Finally, we denote by C the following cone of IR 2 . 

(3.7) C = {t(a, a/K) + (1 - t)(a, Ka):t£ (0, 1), a £ 

Lemma 3.5. Every A £lA is the barycenter of a sequence of laminates of finite order v n £ C 
such that supp v n £lA and 

(3.8) lim / \M\ PK dv n {M) = 00. 

Proof. The proof of Lemma 13.51 follows the strategy of the proof of [5l Proposition 3.10], 
where the particular case of S\ = K and S2 = 1/K is considered. In [5] it is first showed 
that the identity matrix is the barycenter of a sequence of laminates of finite order satisfying 



\Df\ 



VK 



+00 for all balls B £ Q, 
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A 



,1, 




Figure 1. The rank-one connected points A and Q. 



(|3.8p and with support on dC, where C is the cone defined by (|3.7p . The proof is based 
on the construction of the so-called staircase laminates, which was originally made in [8]. 
Then, they extend the result to all other matrices by using the conformal invariance of the 
quasiconvex hull. In our case tl does not enjoy conformal invariance, due to the anisotropy 
of the coefficients <Tj. Therefore, we have to proceed in a different way. 

By slightly modifying the staircase construction in [8], [5] (in fact only a finite number 
of steps at the beginning of the staircase) one can easily show that each point in C can be 
obtained as the barycenter of a sequence of laminates of finite order, satisfying (|3.8p and 
with support on dC. Moreover, by a suitable shift of the support, one can obtain that these 
measures have support in the interior of the cone C. 

Now let A = (ciij) G U. We claim that A is rank-one connected to a diagonal matrix 
Q = (lij) = Qd G C and we conclude the proof. Arguing as in [5, Remark 3.6] it is easy 
to show that Q G U (that is to say, Q belongs to the interior of the quasiconvex hull), and 
that A belongs to a suitable segment [P, Q] still contained in U, i.e., A = tP + (1 — t)Q for 
some r G (0, 1). Since Q G C, Q is the barycenter of a sequence of laminates v n = ^ ^j^A 3 
supported in U and satisfying (I3.8p . The required laminates can then be defined as 



We conclude by proving the claim. From (|3.5I) it follows that ^ £ C. The condition of 
rank-one connectedness reads as 



This is equivalent to the fact that the two rectangles with sides parallel to the axis and 
diagonal Qd-^-d and Q a A a have the same signed area (see Figure [Q) . Notice that the sign of 
the areas is given by the sign of the slope of Qd^-d and of Q a A a , respectively. Define h(Ad, Qd) 
as the signed area of the corresponding rectangle and remark that it is a continuous function. 




(3.9) 



(an 



gii)(«22 - ^22) = a2iai2- 
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Given A, the problem is to find Q d such that 

(3.10) h(A d ,Q d ) = a 2 ia 12 . 

Notice that 

{h(A d ,Q d ),Q d £ C} = [m,oo) 

for a suitable negative m < depending on A d . Therefore, If 021012 > we can always solve 
([3.10p . Assume instead that 021012 < like in Figure [TJ Let h(A d ) be the infimum of h over 
Q d . For a fixed an, it is easy to see that h attains its maximum for 022 = an- In this case, 
the optimal Q d is given by 



Qd = 7^ (an + -j^, Kan + an) , 



and 

a 2 

maxn(^ d ) = ^(1 - 1/K)(K - 1). 

<222 4 

Therefore (|3.10p has a solution whenever 

(3.11) -^i(l-l/^)(^-l)<o 2 io 12 . 
From (13.51) it follows that 



I ct 2 i I < max{Si, S 2 }K5a22, 

and hence 

I012O21I < max{Si, S2}K6a,22\ai2\ < max{Si, S , 2}-fC(5 2 aiia22 < max{5i, S 2 }K8 2 Ka{ 1 . 
By the very definition (|3.4p of 8 we deduce that 

a 2 

-^(1 - 1/K){K - 1) > max{5i, S 2 }K 2 5 2 a 2 u , 
so that (|3.1ip holds, and the proof is completed. 

□ 

4. Two phase Beltrami coefficients 

In the present section we focus on two-phase Beltrami coefficients. In this class, we find 
the ellipticity K mm defined in (|2,25p and we characterize the Beltrami coefficients for which 
K = K mm . From now on, to easy notation, we will omit the dependence on G and H in the 
ellipticity constants. 

4.1. Two-phase Beltrami equation. Let E\ be a measurable subset of Vl and let E 2 '■= 
Q\Ei. Fix {Gi, G2, Hi, H2} C SL sym (2) positive definite (symmetric and with determinant 
one), and consider the functions 

(4.1) G := Xe Gi + X E2 G 2 , H := x E Hi + X E ,H 2 , 

where \ B and Xe 2 are *he characteristic functions of E\ and E 2 , respectively. From (12. lip 
it follows that for G and H of the form (|4.ip , one has 

K = m&x{\gh\ LEi, \gh\ LE 2 } = max{o x /ii , g 2 h 2 }. 

where and hi denote the largest eigenvalue in Ei of G and H, respectively. Set 



K := \/gihig 2 h 2 . 
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Lemma 4.1. The following inequality holds 

K min <k<K. 

Proof. The inequality K < K is trivial. Let us prove that K mm < K. Without loss of 
generality we may assume that g\h\ > g 2 h 2 . Set 



A := J** < 1. 

V 9ihi 

We can have either of the following cases: hi < maxjgi, g 2 , hi, h 2 } or hi = max{gi, g2, hi, h 2 }. 
Suppose we are in the first case. Up to a diagonalization, Gi is of the form 

Gl = ( ? — 

V 91 

We want to use the change of variables (I2.23p . and we recall that g(A, B) and h(A, B) denote 
the maximum eigenvalue of G and H , respectively. We choose 



B 



Vx 

* 



and A = I. Then gi(A,B) = Xgi and g2(A,B) < — Therefore 

A 

gi(A, B)hi (A, B) = K and g 2 (A, B)h 2 {A, B) < K. 

We deduce 

K min <gi(A,B)hi(A,B) = K. 
Suppose now that hi = max{gi, g 2 , hi, h 2 }. Then, after diagonalization of Hi, we choose 

B = I and A = f ^ * ) ' &n< ^ W6 P roceec ^ as before. □ 

Remark 4.2. A direct consequence of Lemma |4. 1 1 is that K mm < K whenever gi hi < g 2 h 2 . 
Proposition 4.3. The following formula for K mm holds: 



(4.2) K mm (G, H) = ^Jg 2 (G^ /2 ,H- 1/2 )h 2 (G^ /2 ,H~ 1/2 ). 

Proof. In view of Lemma 14.14 it is enough to prove that for each A, B € SL(2) we have 

(4.3) g 2 (G~ 1/2 ,H- 1/2 )h 2 (G- 1/2 , H~ 1/2 ) < 9l (A, B)hi(A, B)g 2 {A, B)h 2 (A, B) . 
For this purpose, we show that if Gi = Hi = Id, then for each A, B G SL(2) 

(4.4) g 2 < gi{A,B)g 2 {A,B), 

(4.5) h 2 <hi{A,B)h 2 {A,B). 



Let B e SL{2) and set 
For every v E M? we have 



Gi := B^ B , G 2 := B^ G 2 B . 

|2 ^ 1 ri „. „.\ II d„.I|2 



9i(A,B) 



\v\\ < (Giv,v) = \\Bv\ 
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(4.6) A 1 G\A := ^ _i , A 1 G 2 A : 



and hence 

g 2 {A,B) = sup (G 2 v,v) = sup (G 2 Bv,Bv) > . . , 
IHI<i i|d||<i ffnA^v 

which proves (|4.4|) . The proof of ()4.5p is fully analogous. 

□ 

4.2. Gradient integrability and critical coefficients. In the next proposition we will 
show that if the bound K mm < K is achieved, then Gi and Hi can be simultaneously diago- 
nalized. 

Proposition 4.4. Let G and H be as in (|4.ip and assume that K min = K. Then, there exist 
A,B £ 0(2) such that 

9i \ aT „ a ( i 

52 

(4.7) B^ lB; =(*' °), ^«:=(A° 

Proof. We can always assume that G\ and are as in (|4.6p - (|4.7|) . We prove that, in this 
case, also G 2 is diagonal (For H\ and H 2 we argue exactly in the same way). Set 

B:=G~^, G l :=BGiB = I, G 2 := BG 2 B , 
A := , i^i = AHii = J, # 2 = AH 2 A . 
Since /12 < h\h 2 , g 2 < g±g 2 and recalling Proposition 14.31 we have 

{K min f = gAfok = g 2 h 2 < g 2 h x h 2 < h x h 2 g x g 2 = K 2 , 

where gi and hi are the largest eigenvalues of Gi and Hi. Since K mm = K, all the above 
inequalities are indeed equalities, and in particular g 2 = gig 2 , that implies G 2 diagonal. 

We are left to show that e 2 is the eigenvector associated with g 2 . Arguing by contradiction, 
we assume that 

52 
„ 

Without loss of generality we may suppose that g\ < g 2 and we set 

B := G~^ , Gi := BG X B = I , G 2 := BG 2 B. 

It can be easily checked that gi < gi, that (recall K mm = K) provides the following contra- 
diction 

(K min ) 2 < 9l h ig2 h 2 < 9l h ig2 h 2 = K 2 . 

□ 

We are in a position to show that, for two phase coefficients, Proposition 12.61 is sharp. 

Theorem 4.5. Let G 1} G 2 , H x , H 2 G SL syrn {2), and let K min be as defined in Then 
we have 

i) Let G and H be as in (j4.ip . Every solution f G H\ oc (Q;C) to (]1.4p belongs to 

L Lc(Q c ) f° r ever *y P G %Pk™™); 

ii) There exist G and H as in (14. ip . and a corresponding solution f € H} (0;C) to (ll.4p 
with V/ (5; M 2x2 ) /or every cfofc Bcfi. 



C2 = I n I 
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Proof. The first part of the Theorem is a particular case of Proposition 12. 6| so we pass to the 
proof of ii). By the definition of K mm and by Proposition 14.41 we can always assume that 
G and H are diagonal as in ()4.6p . (|4.7p . with g^hi = K mm . A straightforward computation 
shows that the corresponding cr, defined according to (|2,4p . takes the form 

at := diag [St, K min ) , a 2 := diag (s 2 , , K' 1 <S { <K. 

Therefore, ii) follows from Theorem 13.21 □ 

Remark 4.6. In ii) of Theorem 14.51 we can a ls° enforce that f± satisfies suitable affine 
boundary conditions. 



5. TWO PHASE CONDUCTIVITIES 

In this part we study the gradient summability of cr-harmonic functions corresponding to 
two phase conductivities. Let E\ be a measurable subset of f2 and let E 2 '■= 0\E\. We assume 
that both Ei and E 2 have positive measure. Given positive definite matrices <7i, o~2 € M 2x2 , 
define 

(5-1) cr := x Ex <>\ + X E2 V2- 

Set 

K min = K min^ ._ K min (G (a) , H (a)) . 

5.1. Main results and optimality of the bound (|1.8j) . We can now rephrase Theorem 
14.51 (see also Remark 14.6ft in terms of the coefficient a. 

Theorem 5.1. Let at, o~2 G M 2x2 be positive definite. 

i) Let a be a two phase conductivity as in (15. lft . Every solution u G Hl oc (Q) to (II. ip 
satisfies Vu G L P oc {0) for every p G [2,p K min); 

ii) There exist a as in (|5.ip and (vi,V2) G M 2 such that the solution u G H} (Q) to 

div(aVu) = in 0, 

u{x) = v\X\ + V2X2 on 80, 

satisfies for every ball Bc!J 

\7u\ p x min dx = 00 . 

We are in a position to prove that the bound in (jl.8p is achieved by a suitable conductivity 
a of the type 

(5.2) * = XE 1 (° b b a )+XE 2 (y a b witha = A,6 = ± V / l-A 2 . 

Theorem 5.2. There exist a as in (|5.2p . and a corresponding solution u G H} [Q) of (jl.ip 
with affine boundary conditions such that Vtt ^ L Pk a (U) /or ei>er?/ cfeA; 8cl!, where px x is 
given by (f2T22j) . 
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Proof. By JOJ) we have G, t {a) = I for % = 1, 2, and 

A" 1 vT^A 2 " \ rr / A- 1 -yT^A 2 " 



Therefore if (ct) = K min (cr) = i+v 7 ! 3 ^ = ^ We conclude in view of Theorem O □ 

Finally, we fix the ellipticity A and we characterize the pairs (01,02) corresponding to 
solutions with critical gradient integrability. 

Theorem 5.3. Let px x and p K s v m be as in ()2.22j) . 

i) Let G M (A, fi) be a two phase conductivity as in (|5.ip such that there exists a 
solution u G H 1 ^) of (TT) with Vu $ L^ A (ft;M 2 ); th en takes the following form 

a b \ ( a —b 



= XEx ( _ b a j + XE 2 I b a ), with a = A, 6 = ±y/l- A 2 . 

ii) Let G M S ym (A,fi) 6e a too phase conductivity as in (|5.1|) . suc/i i/iai i/iere exists a 

solution u G fl^fi) 0/ ([Til) Vu ^ L /oc A (S1;M 2 ); i/ien, up to a rotation, tafces 
f/ie following form 

= XBidiag(5'i, A -1 ) + X£ 2 diag(S2, A) , with X < S\, S2 < A -1 . 

Proof i) From Proposition 12.21 it follows that K < 1+V ^~ A2 = K\. On the other hand, i) 
of Theorem 15.11 yields K mm > K\. Lemma 14.11 implies K min = K = K\, thus yielding 
gihi = K mm in both phases. Now apply Proposition 12.31 to conclude that i) holds true, ii) 
Again from Proposition 12.21 Theorem 15.11 and Lemma 14.11 we deduce that K mm = K = \ . 
The thesis follows from Proposition 14.41 □ 

5.2. The explicit formula for K mm . Here we give a direct formula for K min depending 
on 01 and 02. 

Proposition 5.4. Leto~i, 02 G M 2x2 be positive definite. Denote byY>\ andT,2 the symmetric 
part of 01 and 02 respectively, and by d\ and cfo their determinant, 

Sj := af , di := det Sj , z = 1,2. 

Then, 

m + \Jvii l — 4d\d2 n + Vn 2 — 4 



2ja\d 2 



(5.3) K mm -- 

where 

m := (02)11(01)22 + (01)11(02)22 - -((02)12 + (02)21) ((01)12 + (01)21 



n :- 



det 01 + det 2 - - ((01)21 - (01)12) ((02)21 - (02)12) 



If in addition o~±, 02 G M 2 *^, then (|5.3p reduces itself to 



K min = max<; A /-3-, VA2 



— 1/2 — 1/2 

where \\ < A2 are the eigenvalues of a 1 o~20~ 1 
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Proof. From Proposition 14.31 it follows that K mm = \Z*fe^2 where g<i and /12 are the maximum 
eigenvalues of G2 ■= G l G^G^ 1 ^ 2 and H2 := H 1 1 ^ 2 H2H l l ^ 2 respectively. Since by (|2.3I) . 

Gi = 



AdjSj, one has 



Go 



^J£- 1/2 £ 2 £- 1/2 J T . 



171 



The eigenvalues of G2 are solutions to the following equation in A 



det 



So — AS 



0. 



trM 
\fd\d2 



A + l = 0. 



(trM) 2 
did 2 



Set M := £ 2 AdjEi. Since 

det (\H-X2 - XEi) = 
the maximum eigenvalue §2 is defined by 

trM 

Vdid 2 
92 = 

A straightforward computation shows that 

trM = (0-2)11(0-1)22 + (0-1)11(02)22 - -((02)12 + (02)21) ((01)12 + (01)21) 
Similarly, one finds that /i2 is the largest root of the equation 

det(# 2 - XHi) = 0. 

Therefore 

triV + y/(tiN) 2 - 4 
h 2 = - , 

where N := H2 AdjH\. It is easily checked that 

1 



=: m . 



trN 



det 01 + det 02 



^((01)21 - (01)12) ((02)21 - (02)12) 



=: n . 



Now assume that 01 , 02 are symmetric. By (|2,20p we find 52^2 = max ^j- , A2 ^- , where Ai < A2 
are the eigenvalues of 

1 



o 2 



22 



Since by (j2T3]) . G 



V det Oj 

0"2 



Adjoj, one has 

1 V det o 2 1/2 _i 1/2 tT 
- -In o 2 (7 1 J 



(H 2 )22 Vdet 01 
1 1 



(#2 



A 1/2 _l 1/2^ 

:Adj(o/ o 2 o x ' , 



2J22 Jdet(al /2 a^al /2 



, ,/ 1/2 -1 l/2w 1/2 -1 l/2\-l 

det(o 1 / o 2 0/ )(cr 1 / cr 2 o{ ) 



22 
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The eigenvalues of a<i are those of a 



1/2 



a 2 (J l 



1/2 



as soon as we prove that 



i 



(#2)22 = vWl V V"l /2 ) = / 




This follows from the fact that Hi and H2 are diagonal and therefore 




□ 



Remark 5.5. Keeping the notation of Proposition [531 if °2 £ ^sym are positive definite, 
a straightforward computation shows that 



Quasiconformal mappings appear in many branches of mathematics. Only rather recently 
they have shown their power in the theory of composites. In the composite material literature 
one of the typical goals is to determine the so-called "G-closure of a set of conductivities". 
Roughly speaking this means the following. Assume that two matrices, called the conductivity 
of the "phases" and denoted by o\, 02 G -A4(A, £1) are given. Consider a two phase composites, 
i.e. a conductivity a of the form a = o\XE\ + Xe^2 where E\ and E2 are a pair of disjoint 
measurable sets with E\ U E 2 = £1. The task is to find the set of all possible "effective" 
tensors a* that can be obtained by mixing these two phases while letting E\ and E2 vary in 
all the admissible ways. To make this concept precise, one needs to define an appropriate 
concept which is called //-convergence and was invented by Murat and Tartar. This notion 
was a general framework which was necessary to treat the case non-symmetric conductivity 
a which could not be treated by the G-convergence previously introduced by Spagnolo. In 
both cases one can establish compactness results and a notion of closure. We will continue to 
call it G-closure according to tradition even if, in this particular case, one really needs to use 
the //-convergence because the tensor a is not assumed to be symmetric a priori. We refer 
to the recent book of Tartar |18| and reference therein for an extensive treatment. 

In this context, an extensive use of certain special properties of solutions to (jl.ip and 
therefore to (|2.7p . has been made. For an accurate review, we refer to |16j . see Chapter 4. As 
a particularly interesting case, we consider Milton's work computing the so called G-closure 
of a mixture of two materials with arbitrary volume fractions [14] . In the symmetric case, i.e. 
when both phases have a symmetric conductivity, the G-closure was found in the eighties. 
The result has a long history which is reviewed in a very recent work by Francfort and Murat 
[9]. We refer the reader to the reference therein for more details about the original work. 

Milton studied the general case without assuming symmetry. He proved that one can 
recover the G-closure for this case by first reducing the problem to the study of a two- 
phase composite in which, in addition, each phase is symmetric, |14| and Chapter 4.3 in |16| . 
and then applying the results for the symmetric case. Milton explained how his work was 
generalizing previous work by many authors including Keller, Dykhne, Mendelsohn and that, 
in turn, he was inspired by some work of Francfort and Murat and some unpublished work 
by Tartar now available in |18j . Lemma 20.3: in two dimensions "homogenization commutes 
with certain Moebius transformations". Without entering into too many details, we want 
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to emphasize here that the basic ingredients behind these transformations have an elegant 
geometrical counterpart when expressed in terms of the Beltrami equation. 

When a is two-phase, by (|2,3p . so are the matrices H and G. In particular H = H\X\ + 
H2X2- Consider now the equation (|2,ip and make the afhne change of variable / — > F = Af, 
then F satisfies a new equation in which the matrix H is replaced by H A '■= A T HA/(det A). 
Therefore choosing A = H~^R^ with R2 G 50(2) and such that RJH2R =: D2 is diagonal, 
one has 

(6.1) H A = I X i + D 2X 2 

so that H A is diagonal and thus {H A )\2 is identically zero. This in turn implies, by (|2.4p that 
the corresponding conductivity 

G- 1 + (H A ) 12 J 

GA '- {HAh2 

is symmetric. 

We observe, in passing, that applying the same strategy to the domain of / one can 
independently reduce a two-phase G to the form 

(6.2) G B = I X i + G 2X 2 

with G2 a diagonal matrix by a linear transformation x —> Bx. 

In the work of Milton, the "symmetrization" property for a two-phase composites is ob- 
tained as follows. Let A G [0, 1)] and let a G At (A, Q). Set 

(6.3) A= (yc d)' ^A = {aa + bJ)(cI + dJa)- 1 

and let U a = (u^,u^) be any solution to the equation (12. 7ft i.e. cVuJ. = J T Vit^. 



Proposition 6.1. For any two-phase composites, there exists A as in f|6.3|) such that the 
corresponding is symmetric and moreover for some A' G [0, 1) one has G A4(\',£l). 

To continue the argument Milton needs to prove that the G-closure problem relative to 
is mapped one to one into that relative to a. He uses [2] the commutation of the linear 
fractional transformation a — > with homogenization, see also [18], Lemma 20.3. 
Our perspective is to use the following property. 

Proposition 6.2. For any given A as in (16.3P for which G Ai(X',Cl) for some A' G [0, 1), 
there exists 

(6.4, 

such that any solution U-£ A = (u^ a ,u^ a ) to S^Vu^ = J T Vu^ A takes the form 
(6.5) U Sa = A'U a . 

Proof We need to prove that there exist {a', b' , c', d'} such that 

Z A (a'Vul + b'Vul) = J T (c'Vu l a + d'Vu 2 a ), 
which is equivalent to show that 

(o'Sa - c'J T )Vul + (b'E A - d'J T )Vul = 0. 
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We now use the equation cVuJ. = J T Vu 2 and write the previous equation as 

[a'S A - c'J T + (b'E A - d'J T )Ja]Vul = 0. 

One possible solution (actually the only one) is found if the matrix in square brackets is zero 
i.e. if and only if 

a'Sa - c'J T + (b'E A - d'j T )Ja = O S A (a7 + b'Ja) = c'J T + d'a 

S A = (c'J T + d'a){a'I + b'Ja)- 1 
and the latter is equivalent to make the following choice: 



(6.6) A' 



c d 
—b a 



□ 



Proposition 16.21 is the key property to the commuting rule and it is, indeed, a linear change 
of variables in the target space of the underlying quasiregular mapping U = (u,v), solution 
to (g771) . 

Finally one may wonder whether (|6.6f) can be chosen in such a way to have T, A £ A4(X' , Q) 
for some A' > 0. To check this we first note that 

T, A = (acr + bJ)(d + dJa)' 1 = aq " + bJ AdjfcJ + dJa) = 

c z det a + cr 

r (c7 + (cr J J J ) J) = — — — - (cl + dJ 7 <r J ) = 



c 2 det cr + d 2 c 2 det a + d 2 

aca + 6c J + adaJ T a T + bda T aca + bda T + be J + ad det aJ T 
c 2 det a + d 2 c 2 det a + d 2 



It follows that 



s H A + ac + bd s 



Therefore, recalling (|6.6p . the first necessary condition to (|1.2|) can be expressed as follows 
(6.8) c 2 + d 2 > , ac + bd > <S> ac + bd > <S> det A' > 0. 

Now we need to consider S^ 1 . 

Y,- 1 = <d + dJa)(aa + bJY 1 = c/ + dJa Adj(a<7 + 6 J) = 

a" 1 det a + b i 

cl + dJa t t t xT , t t\ (cl + dJa)(aJa T J T + bJ T ) 
'aJa 1 J 1 + bJ 1 ) - 



a 2 det a + 6 2 a 2 det a + o 2 

acJa T J T + adJaJa T J T + bcJ T + bdJaJ T _ 
a 2 det a + 6 2 

acJa T J T + bdJaJ T + adJ T det cr + 6cJ T 
a 2 det a + 6 2 

It follows that 

(6.9) (S^ 1 ) 5 = J£±™ Jo8j*. 

' K A ' a 2 det cr + 6 2 
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Therefore the second necessary condition to (| 1 .3[) is expressed as follows 

(6.10) a 2 + b 2 >0, ac + bd > <t=^ ac + bd > det/4' > 0. 
Putting (j6.8j) and ()6.10|) together we obtain 

(6.11) £ A eA4(A',Sl) for some A' > det A' > . 

Again, this fact has a clear interpretation in the language of the Beltrami system, recalling 
that A' represents a linear change of variables in the target space and that ellipticity in this 
context is measured according to Proposition 12.11 
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